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We propose a scheme for quantum estimation by means of parametric amplification in circuit 
Quantum Electrodynamics. The modulation of a SQUID interrupting a superconducting waveguide 
transforms an initial thermal two-mode squeezed state in such a way that the new state is sensitive 
to the features of the parametric amplifier. We find the optimal initial parameters which maximize 
the Quantum Fisher Information. In order to achieve a large number of independent measurements 
we propose to use an array of non-interacting resonators. We show that the combination of both 
large QFI and large number of measurements enables, in principle, the use of this setup for Quantum 
Metrology applications. 


1. INTRODUCTION 

Circuit Quantum Electrodynamics (QED) [TJ [5] has 
quickly developed in the last decade and is now one of 
the most promising platforms for quantum technologies 
such as quantum computers [m and quantum simula¬ 
tors [51 IS] due to, among other reasons, the high level 
of controlability and scalability that can be achieved. In 
circuit QED both superconducting qubits and electro¬ 
magnetic radiation can be controlled and manipulated to 
an extent that goes beyond some of the standard restric¬ 
tions in other platforms of quantum optics and quantum 
information. 

An important area in quantum technologies is devoted 
to the emergent field of quantum metrology [7] which 
aims at improving the precision of measurement devices 
by exploiting quantum features such as entanglement and 
squeezing in phase estimation protocols. Applications 
are critical and diverse, ranging from the use of squeezed 
states in gravitational wave detection with laser interfer¬ 
ometers |5] to the notion of a global network of quantum 
clocks |5], among many others. 

In this work we propose a scheme for quantum es¬ 
timation in circuit QED. We consider a superconduct¬ 
ing transmission line interrupted by a superconducting 
quantum interference device (SQUID). This technology 
resembles the one employed in the observation of the Dy¬ 
namical Casimir Effect [10]. However, instead of an ini¬ 
tial vacuum we consider the preparation of more general 
initial states- in particular thermal two-mode squeezed 
states- by means of an additional transmission line. The 
modulation of the SQUID transforms this initial state in 
such a way that it becomes dependent on the parame¬ 
ters of the modulating magnetic field. Thus the param¬ 
eters of the magnetic field can be estimated by means of 
phase estimation techniques. We compute the Quantum 
Fisher Information (QFI) and maximize it over the set 
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of considered states in order to determine the optimal 
initial estate for quantum estimation. In order to max¬ 
imise as well the number of independent measurements 
and accordingly the precision, instead of considering a 
single superconducting resonator we propose the use of a 
large array of non-interacting cavities [5]. We show that 
good precision can be achieved for realistic experimental 
parameters. We discuss possible applications of these re¬ 
sults, which include accurate frequency measurements or 
highly precise measurements of magnetic flux variations 
threading the SQUID. 

The structnre of the paper is the following. In Section 
2 we introduce our model and show how to maximize the 
QFI of the electromagnetic field state confined within a 
single superconducting resonator. In section 3 we show 
how these results enable the use of an array of resonators 
for quantum estimation, discussing some potential appli¬ 
cations. We conclude in Section 4 with a summary of our 
results 


2. QUANTUM FISHER INFORMATION IN A 

SINGLE SUPERCONDUCTING WAVEGUIDE 

We will consider a large array of superconducting res¬ 
onators [5] consisting of superconducting waveguides ter¬ 
minated by SQUIDs [B dH with mutual interactions 
controlled by additional SQUIDs m- As we will see in 
detail below we want to achieve a large number of inde¬ 
pendent measurements of a particular quantum state of 
the electromagnetic field. To this end, the SQUIDs can 
be tuned in a such a way that the resonators are non¬ 
interacting m- Therefore, we can focus in the dynamics 
of a single superconducting waveguide as we will do in 
the following. 

In order to exploit squeezing and entanglement we con¬ 
sider the preparation of a two-mode squeezed state as ini¬ 
tial state. This can be achieved by connecting the trans¬ 
mission line to an auxiliary line terminated by an array of 
three SQUIDs which provide a Kerr medium that can be 
used as a parametric amplifier -this has been used to ex¬ 
perimentally generate two-mode squeezed states within a 
single transmission line M- We consider as well a non- 
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zero small temperature characterized by a small number 
of thermal photons 

We will describe the field dynamics by means of the 
covariance matrix V. Using the same convention as in 
which assumes zero displacement without any loss 
of generality, we have Vap = | {RaRp + RpRa ), where 
R7 = (( 7 _,p_, g+,p+) is a vector with the quadratures as 
elements: q± = ( 6 ±+5^)/-\/2 and = — 1 ( 6 ± — 5^)/-\/2, 

given in terms of the creation 6 ^ and annihilation b± 
operators of the two modes of interest +, —. 

The initial state is then described by the covariance 
matrix of a thermal two-mode squeezed state: 


V 

A 

B 


IfAB 
2 \B A 

cosh(2r)(l -I- 2n**')]l, 
sinh(2r)(l -|- 2n*^) cos(0)(Tz 
sinh(2r)(l -|- 2n*^) sin(6*)tTa; 


where r and 0 define the complex squeezing parameter 
X = re*® and are standard Pauli matrices. 


The aim now is to transform this initial state under the 
parametric amplification process induced by the modu¬ 
lation of the SQUID that terminates the waveguide. If 
we add a weak harmonic drive to the SQUID charac¬ 
terized by a frequency Wd and a normalised amplitude 
e, then the field quadratures are transformed as follows 
[H]: g± = -{qo± + fP0T)> P± = -{Po± + /9ot) where 
the small parameter / is: 


/ 


2v 


( 1 ) 


Left is an effective length that describes the boundary 
conditions that the SQUID provides to the flux field while 
V is the speed of light along the waveguide. We are as¬ 
suming that the frequencies of the modes are very close 

w+ ~ uj- ~ Wd/2. 

Under this transformation and considering only up to 
linear terms in n*** and up to quadratic terms in / we 
obtain the transformed covariance matrix V of the state: 


A = cosh(2r)(l -I- 2n*^)(l -I--I- 2/tanh(2r) sin(d))l 
B = sinh(2r)(l -|- 2n‘’')(l — /^) cos{0)az 
+ [ 2 /cosh( 2 r)(l-|- 2 n*^)-|- 

{I + + 271***) sinh( 2 r) sin( 0 )]cra;. ( 2 ) 

Our main aim is to analyze the sensitivity of the state 
in Eq. ([^ with respect to the parameter /. In order to 
achieve this goal we consider the QFI, which provides a 
bound on the error of the estimation of the parameter. 
Therefore, we will seek to maximize the QFI. 


2.1. Single-mode reduced covariance matrix 

First, let us analyze the case in which we try to esti¬ 
mate the parameter by means of measurements over only 
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FIG. 1: (Color online) Single-mode QFI (Hf) vs. r and 6 
for = 8 • 10“^ and / ~ 0.02. The QFI is maximized at 
0 — 7r/2 and 9 = 3n/2 and grows monotonically with r. 


one mode. Therefore, we consider the reduced single¬ 
mode covariance matrix, that is A. 

The QFI for estimation of a parameter r using a single¬ 
mode Gaussian state a is given in m and for zero dis¬ 
placement reduces to: 


_ l Tr[{a-\T)a'{T)f] ^ P'jrf 

" 2 1 + P(r)2 ^ 1-P(r)4’ 


(3) 


where P = l/(4VDet cr) is the purity of the state and Det 
stands for the determinant of a matrix. In our case a is 
the reduced matrix A. The prime indicates a derivative 
with respect to the parameter r (e.g. P'{t) = drP). For 
our purposes this parameter will be /. 

In Fig. ([^ we plot Hf with respect to r and 9, using 
realistic experimental parameters n*** = 8 • 10“^-which 
corresponds to a temperature T = 50 mK-, = 27r x 
10 GHz, Leff = 0.4 mm and e = 0.25 [TUI [II]. We find 
that the QFI oscillates with 9 in such a way that the 
maximum is reached ai 9 — 7r/2, 37r/2 while the minimum 
is at 0 = 0, TT. The QFI grows significantly with the value 
of the squeezing parameter r -as expected- which in the 
figure is plotted in a realistic range r < 2 m- 

In Fig. ([^ we choose the optimal value 9 = TTj2 and 
plot the single-mode QFI vs. r and /. We see that the 
QFI slightly decreases with /, while the growth of the 
QFI with r is observed at any value of /. We have also 
found that the QFI is highly insensitive to the value of 
within the perturbative regime that we are considering 
here. 
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FIG. 2: (Color online) Single-mode QFI {Hf) vs r and / for 
n*** = 8 • 10“^ and the optimal value 0 = ty 12. The QFI 
slightly diminishes with / and grows dramatically with r. 


FIG. 3: (Color online) Two-mode QFI {Hf) vs. r and 6 for 
n**' = 8 • 10~^ and / ~ 0.02. The optimal values are 0 — Q 
and TT while the QFI grows monotonically with r. The QFI is 
three orders of magnitude larger than the single-mode case. 



2.2. Full two-mode state 

Now we analyze the case in which the full two-mode co- 
variance matrix is used for the estimation protocol. The 
two-mode QFI Ht with respect to a parameter r can 
be computed by means of the Uhlmann fidelity in the 
following way HU: 


Ht = —2 lim 


d^F 


s^o d^S 


(4) 


where F is the Uhlmann fidelity given by [TH]: 

F = l/(-\/r -I- '/a — \J (-v/r -I- a/A)^ — T ) and where 

r = 16Det[Q(Vi/2)Q(V2/2) - ^],A = 16Det[(Vi -h 

iQ)/2]Det[(V2 +tQ)/2],T = Det[(Vi -h V2)/2] Q being 
the symplectic form Q = i<Jx (B i<Jx, and the covariance 
matrices Vi and V 2 only differ in an infinitesimal varia¬ 
tion of the parameter of interest, that is Vi depends on r 
while depends on r-|-(5. Thus, in our case, Vi is given 
by Eq.(0 and V 2 is obtained by replacing / by / -I-15. 

Putting all together we are able to find a simple an¬ 
alytical expression for the leading order in perturbation 
theory of the two-mode QFI: 

F[f = 4[sinh^(2r) cos^(0)(l-I-4/^ — dn***) 
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FIG. 4: (Color online) Two-mode QFI {Hf) vs. r and for 
the optimal value 0 = 0 and / « 0.02. We see that thermal 
noise slightly degrades the QFI. 


than the best scenario in the single-mode case. The QFI 
oscillates with 0, but in this case 0 = 0, tt are the optimal 
values. As expected, the QFI grows drastically with r. In 
Figs. I^andj^we see that the QFI is almost independent 
of the value of / and slightly decreases with n}^ -in both 
cases for the optimal value of 6. 

In summary, the strategy to maximize the QFI would 
be to consider the joint two-mode state, to choose an 
optimal value of 0 = 0, tt and to achieve a squeezing 
parameter r as large as possible given the experimental 
limitations. 


In Figj^we plot this two-mode QFI vs r and 0 for the 
same parameters as in We see that the QFI for the 
optimal parameters is three orders of magnitude larger 


In the next section, we will see how all the above is 
related with the error in the measurement of physical 
magnitudes. 
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vs. r and / assuming M = 10^, showing that E can be 
extremely small for the largest values of r. 
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FIG. 5: (Color online) Two-mode QFI {Hf) vs. r and 
/ for the optimal value 6 = 0 and the experimental one 
n*'' = 8-10“® -corresponding to T = 50 mK.The QFI is highly 
insensitive to the value of / while grows significantly with r. 
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3. QUANTUM ESTIMATION OF PHYSICAL 
PARAMETERS IN CIRCUIT QED ARRAYS 

The relation between the optimal uncertainty in the 
estimation of / and the QFI is governed by the quantum 
Cramer-Rao bound [3: 


A/> 




— 1 = 0.0092 
-1 = 0.0184 

1 = 0.0276 

- 1 = 0.0368 


-1 = 0.0092 

-1 = 0.0184 

1 = 0.0276 
-1 = 0.0368 


where M is the number of independent measurements 
performed on the state. There always exists an optimal 
measurement strategy that saturates the bound. 

In the previous section, we have analyzed how to max¬ 
imize the QFI. In the following we discuss how to maxi¬ 
mize the number of measurements. 

To this end we consider a large array of superconduct¬ 
ing resonators mm in which additional SQUIDs control 
the interaction between each resonator m- In particu¬ 
lar, these SQUIDs can be tuned in order to switch off the 
coupling m and to obtain a lattice of non-interacting 
resonators. Thus we can assume that we have a large 
number of copies of the same individual superconducting 
resonator. In this way, we can use a large number M in 
Eq. ([^ in order to mimimize the error. In particular, 
a number of 10^ resonators seems to be within reach of 
current technology [5]. 

In Fig. 1^ we show the number of measurements re¬ 
quired to achieve a relative error E = Af/f <0.1 for 
different values of the parameters r and /. We see that 
for large enough values of r, M can be comparable to the 
desired reference value 10^, although the number grows 
for the lowest values of /. Indeed, in Fig. [^we plot E 



FIG. 6: (Golor online) Number of independent measurements 
M needed to achieve a relative error E < 0.1 in the estimation 
of / vs. r, for different values of /, = 8 • 10“® and 6 = 0. 


The above results show that we can estimate with 
high precision the value of /, that is the degree of ad- 
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FIG. 7: (Color online) Error (E = A///) for different values 
of r and / assuming M — 10^, n*'' = 8 • 10“^ and 9 = 0. 


ditional squeezing introduced by the parametric ampli¬ 
fication. Perhaps more interestingly from the physical 
viewpoint, notice that / is a product of several physi¬ 
cal magnitudes ojd, v and L^g - Eq. ([^. Thus we can 
relate the uncertainty in the estimation of f with their re¬ 
spective uncertainties via the standard error propagation 
formula: 


/ 


dv f 


dujd f 


AL 


( — )2 + (^)2 + (. 

V ujd L 


Y, 


where in the last line we have already used Eq. ([^. 
Assuming that we can consider the scenario in which 
we have good control over all the variables but one - 
for instance, by means of a careful calibration process 
so that we can account for the uncertainties in the 
other variables as systematic errors- then we have that: 
E = A/// = Ax/x -where we are denoting by x the 
mentioned magnitude of interest. That is, the relative 
error in the estimation of / is equal to the relative error 
in the estimation of x. Since we have already shown that 
we are able to achieve very high precision in the estima¬ 
tion of /, this entails that the same conclusion can be 
extended to the estimation of Wd, v and L^g -provided 
that it is possible to realise an experimental scenario in 
which we have control over all the variables but one. 

Moreover, Lgff can be related to the magnetic flux 
threading the SQUID. Indeed, the effective length is 
given by 


Eeg 


/ 4 ’ 0\2 1 

'eJTo 


(7) 


where is the magnetic flux quantum, Lq the induc¬ 
tance per unit length of the superconducting waveg¬ 
uide and Ej the flux-dependent Josephson energy of the 
SQUID, given by: 




( 8 ) 


where Ic is the SQUID critical current and (j)ext tbe ex¬ 
ternal magnetic flux. By combining Eqs. Q and ([^ and 
using an error propagation formula similar to Eq. we 
find: 


ALeff 

Leg 


= Trtan 


(t>0 


^4^ ext 

4>0 


(9) 


In the experiments, 4>ext — 0.35(^o and thus we hnd that 
ALeg/Lgg ~ 2A(^ext/<(>ext ■ Therefore, a precise estima¬ 
tion of / entails a precise estimation of the magnetic flux. 
Of course, highly accurate and sensitive magnetometers 
are already available and indeed SQUIDs are well-known 
as ultrasensitive magnetometers. A thorough investiga¬ 
tion on whether our quantum metrology techniques can 
actually improve on the current state of the art lies be¬ 
yond the scope of this work. 


4. CONCLUSIONS 

We have shown that quantum metrology tools can be 
used to estimate physical variables in a circuit QED sce¬ 
nario. In particular, we have considered a superconduct¬ 
ing waveguide interrupted by a SQUID, where modula¬ 
tion of the magnetic field which threads the SQUID acts 
as a parametric amplifier. 

We start from a thermal two-mode squeezed state 
where the squeezing is characterized by the parameters r 
and 0. We find the optimal parameters that maximize the 
QFI and therefore the initial state which is more suitable 
for quantum phase estimation. After computing both 
the QFI of the full state and the reduced single-mode 
state, we conclude that the best strategy is to consider 
the full two-mode state where the initial parameters are 
0 = 0, TT and r as large as allowed by the experimen¬ 
tal limitations. In order to achieve a large number of 
independent measurements, we propose to use a large 
array of superconducting waveguides, where additional 
SQUIDs control the interaction strength in order to en¬ 
sure a large number of non-interacting superconducting 
resonators, providing copies of the single-resonator sys¬ 
tem. We show that the combination of large QFI and 
large independent measurements enables a precise esti¬ 
mation of the parameter /, which characterizes the pro¬ 
cess of parametric amplihcation. This can be used for a 
precise estimation of the physical magnitudes involved in 
the dehnition of /, for instance the magnetic flux thread¬ 
ing the SQUID. 
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